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Abstract. We calculate the sensitivity of future cosmic microwave background 
probes and large scale structure measurements from galaxy redshift surveys to the 
neutrino mass. We hnd that, for minimal models with few parameters, a measurement 
of the matter power spectrum over a large range of redshifts has more constraining 
power than a single measurement at low redshifts. However, this improvement in 
sensitivity does not extend to larger models. We also quantify how the non-Gaussian 
nature of the posterior distribution function with respect to the individual cosmological 
parameter influences such quantities as the sensitivity and the detection threshold. For 
realistic assumptions about future large scale structure data, the minimum detectable 
neutrino mass at 95 % C.L. is about 0.05 eV in the context of a minimal 8-parameter 
cosmological model. In a more general model framework, however, the detection 
threshold can increase by as much as a factor of three. 
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1. Introduction 

Neutrinos differ from other types of energy density content in our universe in two 
significant ways: They are essentially non-interacting during the entire epoch of 
structure formation, and transit from a relativistic to a non-relativistic particle species 
also in the same epoch. The first aspect means that the dynamics of neutrinos cannot be 
described by simple fluid equations; effectively, the neutrino fluid has infinite viscosity. 
The second aspect means that neutrinos free-stream a finite distance before becoming 
non-relativistic; this introduces a new scale, the neutrino free-streaming scale. 

These unique features imply that cosmology can be used as a very efficient tool to 
constrain the neutrino mass. For neutrinos of sub-eV masses, the signature of massive 
neutrinos on the large scale structure (LSS) power spectrum, i.e., the gradual damping 
of fluctuation power, falls between wavenumbers k ~ 0.01 h Mpc -1 and 0.5 h Mpc" 1 , a 
range conveniently probed by current cosmological probes. Indeed, shape information 
from presently available measurements of the LSS power spectrum, when combined 
with observations of the cosmic microwave background (CMB) anisotropies, can already 
constrain the upper limit on the sum of all neutrino mass eigenstates m v to lie in the 
range 0.2 — > 1 eV in the context of the concordance ACDM model. The precise value 
of the limit depends strongly on the data sets used, as well as on the number of free 
cosmological parameters employed in the analysis (see [1-10] for a selection of recent 
papers on the subject). 

Furthermore, because the fraction of highly non-relativistic neutrinos changes with 
time, and only highly non-relativistic neutrinos can cluster, the impact of neutrinos on 
structure formation depends not only on scale, but also on time. For neutrinos that are 
non-relativistic today, the LSS power spectrum in terms of wavenumber k and redshift 
z is approximately given by [11] 

^ 2 



P(k,z)=P CDM (k,0) x 



.a + %(0), 

.(l + %(0). 



k < k 



fs 



2-6/„/5 

(W,)%(0)/a nr ]- 6 ^ 5 k> k is 



(1.1) 

where /„ = Q u /Q m , a nr = a(m u = 3T U ) is the scale factor at which neutrinos become 
non-relativistic, and g(z) quantifies the decay rate of the gravitational potential. Save 
for the case of f2 tota i = Vt m — 1, g(z) is a time- dependent function, so that the power 
spectrum suppression due to neutrino free-streaming on scales below the free-streaming 
length Af s = 1-njk^ is also time-dependent. Figure 1 shows the suppression in the case 
of m v = 0.08 eV relative to m v = eV as a function of k and z. On large scales 
the two spectra are identical, while on small scales the suppression matches equation 
(1.1). The rapid oscillations in P(k, z)/P m „ =0 (k, z) can be attributed to baryon acoustic 
oscillations (BAO) in the power spectra. 

Clearly, cosmological probes that tap the power spectrum information in terms of 
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Figure 1. The power spectrum suppression factor P(k, z) / P mv =o(k, z) as a function 
of wavenumber k and redshift z for a neutrino mass of ^TJ m v = 0.08 eV. The shadings 
are from 0.96 (darkest) to 1.00 (lightest) in steps of 0.005. 



both the wavenumber k and the redshift z are a priori extremely sensitive to the neutrino 
mass. One notable example is weak gravitational lensing of distant galaxies, for which 
the impact of massive neutrinos has been studied extensively (see, e.g., [12-15]). 

Another example is high-redshift galaxy surveys. The power of combining CMB 
data with LSS measurements from low-redshift galaxy surveys to probe neutrino masses 
has, since the pioneering work of [16], been discussed in many previous studies (see, e.g., 
[17, 18]). In contrast, the advantage of using combinations of both low- and high-redshift 
surveys has, apart from the study of [19], received relatively little attention so far. 

The purpose of the present paper is to explore this last possibility in detail. A 
number of high-redshift galaxy surveys have been proposed for the next decade and 
beyond. The Wide Field Multiple Object Spectrograph (WFMOS) survey proposed for 
the Subaru 8m telescope, for example, plans to observe at redshifts z ~ 1 and z ~ 3 
[20]. The Hobby- Eberly Telescope Dark Energy Experiment (HETDEX) will similarly 
look out to z ~ 4 [21]. Further down in time, the Cosmic Inflation Probe (CIP) mission 
will observe at even higher redshifts in space [22]. 

For simplicity we will adopt the survey set-ups used in [19]. However, our analysis 
differs from [19] in several important ways: (i) We use a simulation-based method 
for our parameter error forecast, which has been shown to have many advantages 
over the popular Fisher matrix approach [23]. (ii) We derive our constraints within 
more general cosmological frameworks that include uncertainties in, e.g., dark energy 
equation of state, and other parameters degenerate with the neutrino mass. This has 
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important implications for the sensitivity to m u and its detection threshold, (iii) We 
opt to work with spherically-averaged, one-dimensional power spectra, instead of two- 
dimensional spectra which in principle contain additional information from geometrical 
and redshift-space distortions. However, in contrast to the smoothed spectrum analysis 
of [19], we include also information from baryon acoustic oscillations. Combination of 
geometrical/redshift effects and BAO will likely lead to even more powerful parameter 
constraints [24,25]. (iv) Several inconsistent assumptions in [19], particularly their 
neglecting the effects of neutrino masses in the calculation of the CMB anisotropies, are 
rectified. 

Lastly, let us note that, aside from the ability to probe the power spectrum 
evolution, going to high redshifts has the advantage of reducing nonlinearities. At 
low redshifts (z < 0.5), uncertainties in the power spectrum in the nonlinear regime 
require that we restrict the use of galaxy survey data to wavenumbers not exceeding 
k ~ 0.1 h Mpc -1 . Indeed, for neutrino masses below 0.1 eV, the effects of nonlinearity 
are already comparable to the damping of fluctuation power due to neutrino free- 
streaming on scales around 0.05 h Mpc -1 at z — 0. At higher redshifts, however, 
nonlinear effects are less prominent, and it is possible to probe safely the LSS power 
spectrum at much larger values of k. This simple fact has important consequences for 
constraints on neutrino masses. 

The paper is organised as follows. In section 2 we describe the two fiducial models 
to be used in our investigation of the sensitivity of future galaxy redshift surveys. In 
section 3 we describe the error forecast procedure: how mock galaxy power spectra are 
generated, the mock survey parameters we use, and the construction of the likelihood 
function. We give in section 4 the predicted sensitivities to various parameters in both 
cosmo logical models, and we summarise our findings in section 5. 

2. Fiducial models 

We consider two spatially flat cosmological models with eight and ten free parameters 
respectively, 

© 8 = {u b , u c , m u , N v , h, t, A s , n s }, 

©io = {u b ,u c ,^m v ,N v ,h,T,A s ,n s ,a s ,w}, (2.1) 

where u) b = Q b h 2 is the baryon density, u c = Vt c h 2 the cold dark matter density, X m v the 
sum of neutrino masses, N v the effective number of fermionic degrees of freedom during 
radiation domination, h the Hubble parameter, r the optical depth to reionisation, 
A s the scalar perturbation amplitude, n s the scalar spectral index, a s the running of 
the scalar index defined at a pivot scale of k — 0.05 Mpc -1 , and w the dark energy 
equation of state parameter. We consider two different models because nonstandard 
parameters such as a s and w are known to be degenerate with J2 m u (e.g., [26]). It is 
interesting to consider whether future galaxy redshift surveys are actually able to break 
these degeneracies. 
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Table 1. Fiducial values of the cosmological model parameters. 



Model parameter Fiducial value 





0.0223 




0.105 




0.0 




3.04 


h 


0.7 


T 


0.09 


A s 


2.3 x 10~ 9 


n s 


0.95 


a s 


0.0 


w 


-1 


n A 


0.74 



In our definition, N v enters the present-day energy density as 

2 N u m u T,m u 

\l v h = = . 

93 eV 93 eV 

During the radiation-domination epoch the total energy density is given by 
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where T 7 and T u are the photon and neutrino temperatures respectively. Our definition 
of N v differs from that in [19]: they use the definition (2.2) for the present-day neutrino 
energy density, but fix N v at 3.04 in equation (2.3). 

The dark energy density Q\ can be derived from the above parameters using 
f^A — 1 — Q u — — ^c- The fiducial values of our model parameters are given in 
Table 1. 



3. Mock data and parameter error forecast 

We use a simulation-based method to estimate cosmological parameter errors from 
various combinations of future CMB and galaxy redshift surveys. Synthetic data are 
generated according to the experimental specifications of the cosmological probe of 
interest, and then analysed using a Markov Chain Monte Carlo (MCMC) package such as 
CosmoMC [27, 28] . The many advantages of MCMC-based forecasts over the popular 
Fisher matrix approach are discussed in detail in [23]. Here, it suffices to reiterate 
that the MCMC method, which probes the entire likelihood hypersurface, generally 
provides more reliable results than does a Fisher matrix analysis based on estimating 
the likelihood curvature around the best-fit point. 

Reference [23] describes in detail an MCMC parameter forecast using synthetic 
CMB data. We extend the analysis to include also synthetic LSS data from galaxy 
redshift surveys as follows. 
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3.1. Galaxy power spectrum 



Galaxy redshift surveys measure the correlation spectrum of galaxy number density 
fluctuations, P g {k). In turn, this correlation spectrum is related to the underlying 
matter power spectrum P(k) via 

P g (k) = b 2 P(k), (3.1) 

where the bias parameter b varies according to the galaxy type targeted by the survey 
at hand, but is generally expected to be independent of k in the linear regime. 

Given some survey design and restricting the analysis to the linear regime, one 
can expect to measure the galaxy power spectrum P g (k) at a range of wavenumbers 
k € [fcmin, fcmax] up to a statistical uncertainty of [29] 

1 
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(3.3) 



is the effective volume of the survey, with n g {r) the expectation value of the galaxy 
number density at coordinate r. When the condition n g P g (k) 3> 1 holds, V e g is 
equivalent to the actual volume of the survey. 

The quantity A In k is the bin size at k in In A;-space. In general we use a very fine 
binning scheme to capture as much of the features in P g (k) as possible, including wiggles 
from baryon acoustic oscillations (BAO) at k > 0.01 Mpc -1 . The effect of smearing due 
to finite window functions, and hence a partial loss of the BAO features in a realistically 
reconstructed power spectrum, will be discussed in section 3.4. 



3.2. Mock survey parameters 

Following [19], we consider galaxy surveys probing three different ranges of redshifts: 

• Gl: 0.5 < z < 2, ground-based, f2 sky = 1500 deg 2 , 

• G2 : 2 < z < 4, ground-based, f2 sky = 300 deg 2 , and 

• SG : 4 < z < 6, space-based, f2 sky = 300 deg 2 , 

where f2 sky denotes the survey's sky coverage. Note that our Gl survey has five times 
the sky coverage (and hence survey volume) of the survey of the same name in [19]. We 
further subdivide each survey by redshift, although we have also checked that, within 
our one-dimensional (i.e., spherically-averaged spectrum) approach, subdivision or not 
makes essentially no difference to the final results. Table 2 shows the central redshift z c , 
effective volume V^f, and other specifications for each of these mock surveys and their 
subdivisions. As a rough guideline, the specifications of Gl are similar to those of the 
z ~ 1 WFMOS survey [20], while G2 is akin to HETDEX [21] or the z ~ 3 WFMOS 
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Table 2. Mock galaxy survey specifications. From left to right, z c is the central 
redshift of each survey slice, fc m i n (/c max ) the minimum (maximum) wavenumber probed 
by the slice, V c e the slice's effective volume, n g the comoving number density of 
galaxies, and b the bias factor. See [19] for a discussion of the calculation of b. 



Survey 


Zc 


&min 


^max 




fig 


b 






[10- 3 h Mpc" 1 ] 


[h Mpc" 1 ] 


[h- 3 Gpc 3 ] 


[l0- 3 h 3 Mpc~ 3 ] 




Gl 


0.75 


5.32 


0.14 


1.65 


0.5 


1.22 


(0.5 < z < 2) 


1.25 


4.54 


0.19 


2.65 


0.5 


1.47 




1.75 


4.26 


0.25 


3.20 


0.5 


1.75 


G2 


2.25 


7.15 


0.32 


0.68 


0.5 


2.03 


(2 < z < 4) 


2.75 


7.11 


0.41 


0.69 


0.5 


2.32 




3.25 


7.18 


0.52 


0.67 


0.5 


2.62 




3.75 


7.29 


0.64 


0.64 


0.5 


2.92 


SG 


4 


5.82 


0.71 


1.26 


5 


4 


(4 < z < 6) 


5 


6.03 


1.01 


1.13 


5 


5 




G 


6.24 


1.50 


1.02 


5 


5.5 



survey [20]. The space-based survey SG should be a reasonable approximation of CIP 
[22]. See [19] for a more detailed discussion. 

Note that the range of wavenumbers probed varies from survey to survey, and, 
indeed, from redshift slice to redshift slice. At the low end of the spectrum, the 
slice/survey volume defines the minimum wavenumber k min available for observation, 

fc min = 27r/\/ 1 ff /3 , (3.4) 

so that the largest observable perturbation wavelength does not exceed the slice's 
effective length scale. At the other extreme, the maximum wavenumber A; max is chosen 
according to the criterion that the dimensionless power spectrum, 

A '(A0 = (3.5) 
does not exceed unity in the linear theory. Rather than setting common 

^min ctud A? max 

for all surveys/slices of interest, we choose to adhere to these "natural" limits, because 
these limits reflect better the true strengths and drawbacks of each survey set-up. 

3.3. Mock power spectrum data and the likelihood function 

The first step of an MCMC forecast consists of generating a vector of N observed data 
points 

d=(p^(k 1 ), P° hs (k N ) ) T , (3.6) 

and the corresponding error covariance matrix 

iV = diag( [AP^)] 2 , [AP g (k N )} 2 ) , (3.7) 

for a given fiducial model and survey/slice. For simplicity, we take the data points 
P° bs (/cj) to be equal to the fiducial power spectrum P^ d (ki), with errors given by equation 
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(3.2). A more rigorous analysis would require that we draw a random P° hs (ki) from a 
Gaussian centered on Pg d (ki) with variance [AP g (ki)} 2 at every k i: but we choose not 
to pursue this avenue. 

Since both the signal and noise are Gaussian-distributed, we can construct a 
likelihood function in the form 

1 -{b 2 v-d) T N-\b 2 v-d) 



L oc exp 
where the vector v, 



V=(P th (h), P th (k N ) 



(3.8) 



(3.9) 



denotes theoretical predictions of the matter power spectrum P{k). Since the N data 
points are assumed to be uncorrelated, the likelihood function is equivalent to 



X 



-21nL = ]T 



x r R obs (fci) -b 2 P th (k^ 12 



(3.10) 



AP g (k t 

up to a constant offset. 

Observe how the bias parameter b enters the likelihood function (3.8). In general 
we are not interested in the exact value of b; instead, we construct an effective likelihood 
function L e g by marginalising over b 2 , 

L cS oc J alb 2 n(b 2 ) L. (3.11) 

If the prior n(b 2 ) is flat, the marginalisation can be performed analytically to give 

2 In L eff 

/ N- 1 vv T N- v 



d T AT 1 - 



v T N 



-^—jd + ln^N^v), 



(3.12) 



again, up to a constant offset. 

For multiple surveys/slices, the combined likelihood is simply a product of 
the individual effective likelihood functions. The CMB likelihood can be also be 
incorporated in a similar fashion. 



3.4- Window functions 

Our treatment so far supposes that the power spectrum, including the BAO wiggles, is 
well-sampled. In reality, however, the finite widths of the window functions used in the 
reconstruction of the power spectrum from a galaxy survey will smear out the power in 
Fourier space, so that little of the BAO features remain in the reconstructed spectrum; 
A separate analysis in terms of real space two-point correlation must be performed in 
order to extract the BAO peak. 

It is interesting to compare the sensitivities of future galaxy redshift surveys both 
with and without BAO extraction. As said, the former case is already covered by our 
default treatment. For the latter case, we mimic the effect of a finite window function 
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by smoothing the fiducial power spectrum with a normalised top-hat function in log 

space, 

lnP 9 sm (A;) = J dink' W(k,k',R) lnP g {k'), (3.13) 

where 



W(k, k',R) = 




Ink- R/2 < In k' < In k + R/2, 
otherwise. 



We use R = 0.33, roughly similar to the "la" width of the SDSS LRG window functions 
[30]. Mock data points are then drawn from the smoothed spectrum P^ m (/c) with errors 
given by (3.2). 

3.5. Mock CMB data 

Since the LSS power spectrum is generally fairly featureless, parameter degeneracies 
abound so that measurements from galaxy redshift surveys alone do not place very 
stringent limits on cosmological parameters. In order to break these degeneracies, it 
is common to consider power spectrum measurements together with data from CMB 
observations. In the present work, we consider prospective CMB data from the Planck 
satellite [31]. We assume measurements of the auto and cross correlation power spectra 
of the CMB temperature and E-type polarisation (i.e., TT, TE, EE), up to a multipole 
£ = 2250. Experimental characteristics of the Planck satellite, assuming one year of 
observation, can be found in Table 1.1 of the Planck Bluebook [32]. 

4. Results 

Figure 2 shows the projected ID marginalised 95 % confidence regions for various 
parameters within our 8- and 10-parameter models, assuming data from Planck plus 
different combinations of mock galaxy surveys. 

Consider first the 8-parameter model. We note that there is a marked improvement 
in the projected bounds on J2 m u when galaxy survey data from Gl (0.5 < z < 2) and 
SG (4 < z < 6) are combined. This is particularly interesting because the Gl survey 
has more than twice the effective volume of the SG survey. However, since SG is able 
to access a higher A; max and, being a space-based probe, also suffers from less shot noise, 
the Planck+Gl and Planck+SG bounds on J2 m u are almost identical: J2 m u < 0.13 eV 
(95 % C.L.). Furthermore, the fact that the underlying LSS spectrum is both scale- and 
time-dependent (as demonstrated in Figure 1) means that the Gl and SG mock data 
sets contain in principle different parameter degeneracies. This can be seen in Figure 2: 
the SG survey has a better handle on cu c , N u and n s , while Gl is more constraining for h. 
Combination of these data sets, therefore, can help lift these degeneracies (even if only 
partially). Indeed, the combination of Gl and SG yields a bound of J2 m u < 0.08 eV 
(95 % C.L.), which is 60 % better than that from Gl or SG alone. A similar trend can 
also be seen in the allowed regions of n s , further indicating that we are indeed gaining 
better shape information from the combination. 
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Figure 2. Projected ID marginalised 95 % confidence regions for various cosmological 
parameters within the 8- and 10-parameter models, using various combinations of 
Planck and mock galaxy surveys. In each subplot, the data sets under consideration 
are, from top to bottom, Planck+Gl, Planck+G2, Planck+SG, Planck+Gl+G2, 
Planck+Gl+SG, and Planck+Gl+G2+SG. Solid lines denote constraints in the 
context of the 8-parameter model, while dotted lines refer to the 10-parameter model. 



This bound of Yl, m v < 0.08 eV from Planck+Gl+SG is particularly noteworthy 
if we compare it with what can be achieved in a low redshift survey with an effective 
volume equal to the combined volume of Gl and SG. As an example, we might consider 
the survey "Glxl.45", which has 1.45 times the effective volume of Gl. From volume 
considerations alone, the improvement in the parameter errors from Glxl.45 relative 
to Gl is at most a factor of \/TA5 ~ 1.2. But this is a very optimistic figure because 
the parameter constraints do not depend only on data from galaxy surveys; CMB data, 
too, play a crucial role, and there is no corresponding increase in the CMB sky coverage. 
Using an MCMC analysis, we find that Planck+Gl x 1.45 yields a 95 % upper bound on 
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Figure 3. Same as Figure 2, but with power spectrum smoothing as described in 
section 3.4, i.e., no BAO extraction. 



J2 m u of 0.12 eV. Compared with 0.13 eV from Planck+Gl, the actual improvement is 
negligible and cannot compete with the gain from going to higher redshift s. 

For the 10-parameter model, however, the projected sensitivity to m v is much 
poorer and there is little gain even when both the G2 and SG data sets are added to 
the analysis. In this model the effect of the neutrino mass can be mimicked by changes 
in a combination of other parameters, in this case mainly the dark matter density u c , 
the effective number of neutrino species N u , the dark energy equation of state w, and 
the running of the spectral index a s . The larger k range probed by SG does lead to 
a better constraint on a s . Unfortunately this improvement does not translate to the 
neutrino mass bound. 

We caution however that our treatment here, namely, the use of spherically- 
averaged power spectra instead of two-dimensional spectra with additional geometri- 
cal/redshift effects, may not be making full use of the capacity of high- redshift galaxy 
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Table 3. A summary of Figures 2 and 3: projected ID marginalised 90 % (95 %) 



upper bounds on 


Y] m v , assuming 


£m£ d = (i.e. 


, the "sensit 


ivity"). 


Survey 


8-parameter model 


10-parameter model 


BAO 


No BAO 


BAO 


No BAO 


Planck 


0.62 (0.73) 


0.78 (0.96) 


+Glx5 


0.11 (0.13) 


0.12 (0.14) 


0.17 (0.19) 


0.18 (0.20) 


+G2 


0.15 (0.18) 


0.20 (0.24) 


0.26 (0.30) 


0.30 (0.34) 


+SG 


0.12 (0.14) 


0.17 (0.19) 


0.30 (0.34) 


0.34 (0.38) 


+Glx5+G2 


0.10 (0.12) 


0.11 (0.13) 


0.17 (0.19) 


0.18 (0.21) 


+Glx5+SG 


0.069 (0.080) 


0.086 (0.097) 


0.16 (0.18) 


0.18 (0.21) 


+Glx5+G2+SG 


0.069 (0.080) 


0.087 (0.10) 


0.16 (0.19) 


0.18 (0.20) 



surveys. Previous studies show that the combination of geometrical effects and BAO 
can be a very powerful tool to probe separately the Hubble expansion H(z) and the 
angular diameter distance D^(z) [24,25]. This may be useful for breaking the small 
degeneracy that exists between £ m v and w. 

Figure 3 shows the projected ID marginalised 95 % confidence regions for the same 
8- and 10-parameter models, but now using the smoothed spectra described in section 
3.4. Compared with Figure 2, we see that the inclusion of BAO information does lead 
to better constraints on parameters in general. However, it does not play a role in 
improving in the neutrino mass bound between Planck+Gl+SG and Planck+Gl in 
the 8-parameter model; the improvement seen in Figure 2 remains intact in Figure 3 
even without BAO extraction. Clearly, the improved constraints come from tracing the 
evolution of the LSS spectral shape. 

Lastly, we note the 95 % allowed ranges are in general not centred on the fiducial 
model because of the non-Gaussian nature of the marginalised posterior distribution. 
The asymmetric error bars highlight the importance of using MCMC-based forecasts. 

4-1. Sensitivity 

It is interesting to compare the potential of galaxy redshift surveys to constrain neutrino 
mass with the expected performance of laboratory experiments such as KATRIN. The 
sensitivity of KATRIN, defined as the 90 % upper limit assuming a fiducial neutrino 
mass of zero, is estimated to be 0.2 eV [33]. Derived under the same assumption of 
£ ml d = 0, Table 3 summarises the projected 90 % and 95 % sensitivities for various 
combinations of Planck and mock galaxy surveys, with and without BAO extraction, 
for both the 8- and 10-parameter models. 

4-2. Detection threshold 

In addition to the sensitivity, KATRIN quotes a detection threshold, defined as 
the fiducial neutrino mass at which an Na detection is possible. With its present 
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Table 4. Projected ID marginalised 95 % lower and upper bounds on Y m v using 
Planck+Gl+SG with BAO extraction for various fiducial Y m v values. 



8 parameters 10 parameters 



£m« d [eV] 


Lower 


Upper 


Lower 


Upper 


0.0 


0.0 


0.080 


0.0 


0.19 


0.050 


0.0 


0.11 


0.0 


0.18 


0.093 


0.027 


0.14 


0.0 


0.21 


0.16 


0.11 


0.20 


0.0 


0.26 


0.19 


0.14 


0.22 


0.10 


0.28 


0.22 


0.18 


0.26 


0.14 


0.31 


0.25 


0.21 


0.28 


0.17 


0.33 


0.28 


0.23 


0.31 


0.21 


0.36 


0.37 


0.32 


0.40 


0.30 


0.44 



configuration, KATRIN will achieve a 3a detection of an effective electron neutrino 
mass of 0.3 eV and 5<r for 0.35 eV [33]. 

Calculating such a threshold for cosmological data is more involved because it 
necessitates a scan in fiducial model space. Table 4 and Figure 4 show the results of 
such a scan; we derive the ID marginalised 95 % confidence regions using different values 
of Y rn^f. The projected 95 % detection threshold is then defined as the value of Y 
at which the 95 % lower bound differs from zero.f In the 8-parameter case the threshold 
thus defined lies between 0.05 eV and 0.09 eV, whereas in the the 10-parameter model 
it falls in the range 0.16 eV to 0.19 eV because of the degeneracies between Y m u, ®s 
and w. Comparing these numbers with their corresponding 95 % sensitivities derived 
in the last section (0.08 eV and 0.19 eV respectively), we see that they are very similar, 
as one would intuitively expect. 

Let us now consider what might be expected from a forecast based on the Fisher 
matrix, which assumes Gaussianity in the posterior distribution with respect to the 
model parameters. There are two approaches — and two pitfalls. 

The first approach is to derive a as function of Y m u, ignoring the fact that 
the posterior distribution may be non-Gaussian for certain values of Y m u- Besides 
numerical instabilities (see [23] for a discussion), the a values derived in this manner, and 
hence any confidence regions constructed therefrom, have no meaning in the Bayesian 
context, simply because the posterior distribution is non-Gaussian. 

The second approach is to avoid the non-Gaussian region by evaluating the Fisher 
matrix at a larger fiducial value of Y m u at which the distribution is two-tailed and 
(hopefully) approximately Gaussian. The result is then extrapolated to other regions of 
parameter space. For example, in the case of the 10-parameter model, one might want 

| One-tail limits are calculated when the marginalised posterior P drops below 15 % of its maximum 
at Y m v = 0. This corresponds to — 21nP/P max ~ 4, i.e., roughly 2a in the case of a Gaussian 
distribution. 
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0.50 F 




Fiducial T,m u (eV) 

Figure 4. The projected ID marginalised 95 % confidence range for ^ m„ as a 
function of the fiducial m u for Planck+Gl+G2+SG with BAO extraction. The 
light shaded band denotes the 10-parameter model, while the dark band corresponds 
to the 8-parameter model. 

to evaluate the Fisher matrix at, say, J2 m u = 0.25 eV. Granting that we can overcome 
the issue of numerical instability, we should find a "2a" of approximately 0.08 eV at 
J2 m u d — 0.25 eV (see Table 4). However, the curvature of the likelihood function is 
much larger for large values of the fiducial neutrino mass, so that any such estimate 
of the detection threshold will likely be much too optimistic. Indeed, extrapolating 
2cr = 0.08 eV to lower values of Z)m^ d , one would be led to conclude that the detection 
threshold is 0.08 eV, which is clearly much lower than (0.16 — > 0.19) eV! 

This highlights the importance of using MCMC techniques on mock data instead of 
the Fisher matrix approach in order to get reliable estimates of the potential of future 
experiments. 

5. Conclusions 

We have calculated the sensitivity of future CMB probes and LSS measurements from 
galaxy redshift surveys to the neutrino mass. In particular, we have studied how the 
measurement of the LSS power spectrum at different redshifts can help to constrain 
J2 vn v more efficiently than a single, more precise measurement at low redshifts in some 
minimal cosmological models. 

Along the same line we also comment on the difference between the sensitivity, 
the formal upper limit given an underlying neutrino mass of zero, and the detection 
threshold, the minimum m u measurable at a given confidence level. We have 
quantified the difference between the sensitivity and the detection threshold in specific 
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cases of CMB+LSS measurements. We suggest that this method is the most reliable 
for making error forecasts in cosmology. Furthermore, it has the advantage that the 
results can be directly compared with projected sensitivities and detection thresholds 
of laboratory experiments such as KATRIN that probe directly the absolute neutrino 
mass scale. 

Within the context of a minimal 8-parameter cosmological model, we find that a 
minimum neutrino mass of order 0.05 — > 0.09 eV can be detected at 95 % confidence 
level using Planck and a combination of low and high redshift galaxy surveys, while the 
95 % sensitivity is 0.08 eV. The latest neutrino oscillation data prefer the mass splittings 
Am 2 12 ~ 8 x 10~ 5 eV 2 and Am| ~ 2.5 x 10~ 3 eV 2 [34, 35]. If the lightest eigenstate has 
zero mass, these mass splittings imply J2 m u — 0.06 eV for the normal hierarchy and 
0.1 eV for the inverted hierarchy. With the survey configurations chosen for our study, 
it will not be possible to differentiate definitively between the two hierarchy schemes. 

In more complex cosmological models with more free parameters, the detection 
threshold and sensitivity become even worse. In a 10-parameter model including also 
a running scalar spectral index and a nontrivial equation of state for the dark energy 
as free parameters, Table 4 shows that the sensitivity to m u worsens by a more than 
a factor of two compared to the 8-parameter case. The corresponding increase in the 
detection threshold is threefold so that a positive detection will only be possible if 
neutrino masses are in the quasi-degenerate regime. 

We note that the survey configurations employed in our analysis are not necessarily 
overly optimistic. A very recent study of neutrino mass measurements from Planck 
plus a future wide-field galaxy survey out to z ~ 2 [36] assumes a much larger effective 
survey volume than our Gl survey (roughly comparable to Glxl2). We have checked 
explicitly that the 95 % sensitivity to J2 m u in this case is 0.06 eV for the 8-parameter 
model, a figure that can almost be matched by our 0.08 eV from a combination of low- 
and high-redshift measurements with a much smaller total survey volume. 

Finally we stress that although future high-redshift galaxy surveys will be extremely 
useful for probing small neutrino masses, they are by no means the only source of 
information about the LSS power spectrum. Weak gravitational lensing of distant 
galaxies, for example, will provide a complementary and possibly even more sensitive 
probe of neutrino properties [12-15]. In the more distant future, surveys of the 21- 
cm line from beyond the epoch of reionisation will provide access to the matter power 
spectrum at even higher redshifts. 
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